A noodle, hockey stick, and 
spaghetti plate: a perspective on 


high-resolution paleoclimatology 


David Frank, '* Jan Esper,? Eduardo Zorita? and Rob Wilson“ 


The high-resolution reconstruction of hemispheric-scale temperature variation 
over the past-millennium benchmarks recent warming against more naturally 
driven climate episodes, such as the Little Ice Age and the Medieval Warm 
Period, thereby allowing assessment of the relative efficacies of natural and 
anthropogenic forcing factors. Icons of past temperature variability, as featured 
in the Intergovernmental Panel on Climate Change (IPCC) reports over nearly 
two decades, have changed from a schematic sketch in 1990, to a seemingly 
well-solved story in 2001, to more explicit recognition of significant uncertainties 
in 2007. In this article, we detail the beginning of the movement to reconstruct 
large-scale temperatures, highlight major steps forward, and present our views 
on what remains to be accomplished. Despite significant efforts and progress, 
the spatial representation of reconstructions is limited, and the interannual 
and centennial variation are poorly quantified. Research priorities to reduce 
reconstruction uncertainties and improve future projections, include (1) increasing 
the role of expert assessment in selecting and incorporating the highest quality 
proxy data in reconstructions (2) employing reconstruction ensemble methodology, 
and (3) further improvements of forcing series. We suggest that much of the 
sensitivity in the reconstructions, a topic that has dominated scientific debates, 
can be traced back to the input data. It is perhaps advisable to use fewer, 
but expert-assessed proxy records to reduce errors in future reconstruction 


efforts. © 2010 John Wiley & Sons, Ltd. WIREs Clim Change 2010 1 507-516 


Goa the climate sensitivity! and the 
reconstruction of large-scale temperature? may 
be regarded as key goals to help predict consequences 
of anthropogenic activity. The former quantifies how 
global temperature will change due to shifts in 
external (e.g., solar, volcanic, and anthropogenically 
emitted greenhouse gases and aerosols) forcing. The 
latter characterizes variability of the climate system 
due to changes in natural external forcing and 
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internal climate variability (e.g., ENSO, meridional 
overturning circulation) in absence of significant 
anthropogenic forcing. Both of these objectives serve 
in parallel to quantify climate system behavior. 
A large natural variability may imply a large 
climate sensitivity and underscore the consequences 
of anthropogenic emissions.? However, respecting 
uncertainties in both the temporal course and 
efficacy of solar and volcanic activity, a large 
amplitude over the past millennium may equally 
well indicate a greater role of natural factors in 
dictating past (and future) climate variation.’ In 
addition to the direct assessment of the climate 
sensitivity to a doubling of CO2, the amplitude of past 
temperature change is also intimately linked with the 
modulation of climate drivers of many earth system 
processes, such as the positive feedback between 
temperature and atmospheric CO2 concentrations 
as mediated by a variety of oceanic and terrestrial 
processes.”*° 
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To help address these issues, over the last 
30 years, more than two dozen large-scale reconstruc- 
tions of past temperature have been produced—most 
since the late 1990s. These reconstructions have placed 
the twentieth century warmth within the context of 
preceding climatic phases, such as the Medieval Warm 
Period (MWP) and Little Ice Age (LIA), and arrived at 
the broad consensus of the unusual nature of late 
twentieth century warming. The sheer number of 
reconstructions testifies to their scientific and soci- 
etal importance. But simultaneously, this number of 
reconstructions, and continued efforts for still more, 
clearly demonstrate that not all questions of large- 
scale climate variability are solved. In this article, we 
(1) briefly review the beginning of this reconstruction 
movement, (2) describe the ‘consensus view’ for the 
large-scale temperature evolution over the past mil- 
lennium as illustrated in IPCC reports, and (3) present 
what we have learned in the past decade or so and 
what, in our opinion, still remains to be answered. 
We conclude with our view on research agendas that 
should help reduce remaining uncertainties. 


EARLY RECONSTRUCTIONS 


In a 1979 article, Brian Groveman (doctoral stu- 
dent) and Helmut Landsberg (advisor) presented the 
first quantitative reconstruction of Northern Hemi- 
sphere (NH) annual mean temperatures in Geo- 
physical Research Letters.” The authors compiled 
20 time-series of mostly instrumental and a few proxy 
records (including two tree-ring series) that individ- 
ually correlated between 0.28 and 0.62 with NH 
temperatures, and using a nested multiple regression 
to take into account the different temporal coverage of 
the predictors, provided a reconstruction and associ- 
ated measures of uncertainty back to 1579. Groveman 
and Landsberg called this a ‘short cut method to esti- 
mate NH temperature’ based upon a demonstration 
that nine instrumental stations could sufficiently rep- 
resent the mean of an extensive gridded network of 
instrumental series. Aside from the fact that this short 
publication (2 tables, 1 figure, 10 references) is the first 
quantitative reconstruction of NH temperatures, novel 
methodological aspects of the paper appear to have 
only recently been re-discovered. Such aspects include: 
separate calibration for all predictor subsets, time- 
varying error bars, and mention of the (small) large- 
scale amplitude of reconstructed temperatures derived 
from their approach. This groundbreaking reconstruc- 
tion showed evidence that the Maunder Minimum was 
not appreciably cooler than neighboring periods and 
that the twentieth century was rather warm. 
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Remaining untouched for about a decade, 
Jacoby and D’Arrigo, re-approached quantitative 
reconstructions for the NH and developed the first NH 
reconstruction’ based entirely upon noninstrumental 
sources: tree-rings. Building upon the quantitative tra- 
ditions of dendroclimatology,'® Jacoby and D’Arrigo 
presented calibration/verification statistics and miti- 
gated potential problems due to multicollinearity of 
the predictors by performing principal components 
regression. Even though the suite of techniques to 
preserve long-term climate signals in tree-ring records 
were not so well developed at that time, they pointed 
out that low-frequency variance must be preserved 
in the proxy records for them to contribute usefully 
to global change studies. This reconstruction also 
showed that twentieth century warmth was anoma- 
lous back to at least 1671. 


CLIMATE CHANGE ICONS 


Following in the footsteps of these initial efforts”? 
and bolstered by the increased recognition of the 
relevance of global change research and development 
of new proxy records and techniques, numerous large- 
scale temperature reconstructions in the following two 
decades were produced.!!~*° It was not until the late 
1990s that such reconstructions extended far enough 
back to place recent temperature change in the context 
of the perhaps temporally nearest natural warm 
analog: the MWP. Such analyses are of relevance as 
they allow the forcing and variability responsible for 
natural warm and cool conditions to be understood 
and contrasted with the current and more diverse 
forcing cocktail. However, the IPCC reports took an 
early stance on this topic in the first report published 
in 1990, and followed this up in the 2001 third?® 
and 2007*? fourth assessments. In Figure 1, we show 
snapshots of the evolving ‘consensus view’ of large- 
scale temperature change over the past millennium. 
It should be noted that this figure (and related text) 
only represents a small percentage of the knowledge in 
the various IPCC reports related to past-millennium 
temperature change. For example, the 1995 IPCC 
report displayed a submillennial length record that 
extended back to a.D. 14001? and the 2001 report 
also depicted a harbinger of the 2007 spaghetti plate 
by showing three reconstructions. 1416-18 Past regional 
warming, cooling, and the general level of climatic 
hetero/homogeneity as evidenced by the text of the 
IPCC volumes is also not captured by such figures. 


IPCC 1990: ‘Eurocentric Cartoon’ 


In the first IPCC report, presentation of temperatures 
over the past millennium appears to be a hand-drawn 
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FIGURE 1 | IPCC icons of temperatures over the past A oe 
millennium. Sequence of the pre-industrial to industrial = 
temperatures as expressed in the 1990 (upper), 2001 $ 
(middle), and 2007 (lower) IPCC reports. Dashed lines b 
represent mean temperatures at 1900 in the upper panel and = 
for 1961-1990 in the middle and lower panels. The upper = 
panel was graphically recreated, whereas the middle and (®) 
lower panels are based upon data obtained at the NCDC E 


webpage (http://www.ncdc.noaa.gov/paleo/recons.html). 
Curves illustrate the evolution of the ‘consensus views’ for 


IPCC 1990 


Little ice age 


Medieval warm 
period 


IPCC 2001 


IPCC 2007 


large-scale temperature change in the IPCC reports. See text 
for details. 


sketch of how climate varied. A prominent MWP, 
followed by a cold period—the LIA—are observed 
with clear indications for high natural variability of 
the climate system. It was demonstrated*® that this 
figure was derived from the groundbreaking work 
of H. H. Lamb*! representing temperature change 
over Central England, but the chosen smoothing and 
figure description suggested more of an illustration, 
rather than quantification, of past temperatures. This 
curve presenting a Eurocentric schematic view of 
climate change, left room open for a completely new 
temperature history. 


IPCC 2001: ‘Anthropogenic Admission’ 

The 1990 sketch was replaced 11 years later by a 
quantitative reconstruction for the entire NH by 
Mann et al.'° Annual precision and y-axis in degrees 
Celsius contrasted strongly with the 1990 figure and 
suggested a marked increase of knowledge. The little 
varying, yet slowly declining temperatures over about 
900 years of natural variability, reversed by the strong 
upward trend over the past century, served as a 
symbol for human impact on temperature change. 
The existence of an MWP and LIA were hinted at, but 
only as part of an almost negligible pre-anthropogenic 
temperature change. Even though this reconstruction 
was presented with calibration uncertainties, their 
common omission in subsequent discussion, may have 
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led to an overly optimistic assessment of the degree to 
which past temperatures were understood. 


IPCC 2007: ‘Stepping Back’ 


In the most recent IPCC report, further consideration 
of existing and the development of new reconstruc- 
tions, methodological disputes, and analysis called 
for a retreat from the 2001 position that recon- 
structed temperatures were well understood. Numer- 
ous, smoothed reconstructions in the 2007 report 
testify to significant remaining uncertainty. The super- 
position of the instrumental data still demonstrate 
the extreme warmth of the most recent decades, but 
the wide-scatter in the reconstructions in defining 
pre-industrial temperature variability showed that an 
understanding of natural variability remains poorly 
constrained. The pre-1000 A.D. period hints at a rise of 
temperatures into the MWP, with the convergence of 
most reconstructions toward peak values ~1000 years 
ago suggesting for some sort of large-scale MWP. 
However, the high dispersion of the records away 
from the calibration period mean casts doubt on the 
absolute amplitude of past temperature change. 


WHERE WE STAND TODAY 


The so-called ‘Hockey Stick’? debate*?*> is/was 
centered on the not-only academically relevant task of 
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defining the true course of temperatures over the past 
millennium and trying to resolve discrepancies among 
reconstructions due to data and methodological 
details. Similar scientific debates occurred earlier. 
It was pointed out that the 1979 Groveman and 
Landsberg reconstruction received criticism for sparse 
data representation during the early reconstructed 
period and also for the warmth reconstructed during 
the late seventeenth and early eighteenth centuries.” 
The former criticism could be (and perhaps has been) 
invoked for nearly every reconstruction published to 
date. Without wishing to rehash all such details, in 
the following text we highlight where we broadly find 
that progress has been achieved and where more is 
necessary. 

There is now broad recognition for the impor- 
tance of preserving all possible climatic frequencies in 
proxy data and reconstructions. This is a step forward 
to simply using data that have likely not preserved low- 
frequency signals. Advances include the widespread 
use of Regional Curve Standardized tree-ring data 
to preserve the long-term signals,** and efforts to 
derive the highest possible resolution from records 
typically analyzed at decadal or lower resolution.* 
The most recent generation of reconstructions now 
typically employ only highly resolved data where 
the long-terms trends have been preserved (see, for 
example, Refs 22 and 28), and proxies where the 
low-frequency variability has not been preserved are 
increasingly recognized as having little room to con- 
tribute to understanding long-term change. It should 
be noted that the fact that a proxy record is annually 
resolved and thus allows high-frequency information 
to be discerned, should not automatically disqual- 
ify this record to contribute to the quantification 
of low-frequency variability.2° Such selection would 
unnecessarily result in a significant reduction of pre- 
dictors that have skill across all frequency domains 
and possibly overweight low-resolution series for esti- 
mating long-term trends.7% 4°47 More research is still 
needed to consider the frequency dependent charac- 
teristics of proxies’ signals and achieve appropriate 
high versus low-frequency weighting.4”~°° 

The climatic signal in proxy records is now 
routinely calibrated—and often verified—via compar- 
isons with instrumental records allowing assessment 
of proxy trustworthiness in the recent period. 
Interdisciplinary advances have led and are leading 
to increases in the number of available high-quality 
proxy archives. This includes temporal extensions 
both toward present (allowing better calibration) 
and back in time (allowing longer term perspectives 
on natural climate variation). However, it is also 
recognized that reconstructed values are highly 
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dependent upon the calibration models applied*!*? 
and the instrumental data utilized.°3°4 Repetitive 
testing of proxy data against various meteorological 
seasons and parameters until acceptable correlations 
are found likely results in an optimistic view of proxy 
quality.°> In addition, data overlap and data sparse 
regions in large-scale reconstructions, challenges 
independent and truly large-scale assessment of 
climate variation.** Very few temperature sensitive 
reconstructions are available for the mid-to-low 
latitudes, although on going activities are filling these 
gaps.°°>’ Even the most well-replicated large-scale 
reconstructions tend to have only a handful of proxies 
prior to ~1400.78 Techniques that use co-variance 
information among proxy records and instrumental 
fields attempt to fill some of these gaps, but it often 
is not clear how individual proxies are weighted 
in the final reconstruction.°* Although much recent 
literature is devoted to methodological disputes 
and improvements, sparse and noisy data are likely 
the underlying cause for the high methodological 
sensitivity in the first place. 

The regional to large-scale amplitude remains 
poorly characterized. The wide range of reconstructed 
amplitudes presented in the 2007 IPCC report that 
results from variable preservation of low-frequency 
signals, calibration methodologies, and proxy and 
instrumental noise has not yet been fully reconciled.’ 
New ensemble reconstruction estimates that consider 
and mitigate many of these uncertainties place the 
amplitude for annual NH land and sea-surface temper- 
atures at ~0.7°C between the warmest (1971-2000) 
and coldest (1601-1630) climatological periods of 
the past millennium.” Potential biases of long-term 
instrumental temperature trends for both land**°? 
and sea-surface® temperatures affect the calibra- 
tion of low-frequency signals and the temperature 
amplitude. Challenges in defining appropriate spatial 
and seasonal representativity also result in amplitude 
uncertainties. Apart from the amplitude differences, 
the various large-scale reconstructions still display 
large discrepancies in their time evolution, particularly 
at the interannual and centennial time-scales. 

To varying degrees, all three IPCC icons con- 
form to the broad notion of the MWP, LIA, and recent 
warmth. The horizontal lines provided in the two most 
recent IPCC reports (see Figure 1) readily allow com- 
parison and conclusions for the unprecedented nature 
of current (1961-1990) warmth. Doubt is rather most 
strongly cast upon the magnitude of LIA cold. The 
prominent MWP of Lamb (1965) appears to have 
fueled debate on the spatial extent of this era.?8-32,61-63 
For example, in the 2007 IPCC report eight regional 
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proxy series were purported to illustrate ‘the hetero- 
geneous nature of climate during the “MWP””. How- 
ever, a recent analysis showed that based upon these 
data the IPCC conclusions were nontenable and that 
the MWP was no more heterogeneous than the LIA or 
even the present.®? Rather, it was suggested that too 
few (and noisy) proxy series currently exist to assess 
possible heterogeneity. A recent reconstruction of the 
North Atlantic Oscillation (NAO; see Ref 64) showed 
a persistently positive low-frequency state between 
about 1050 and 1400 and thus helped provide a 
dynamical explanation for European Warmth during 
Medieval times.*° However, globally distributed 
proxy data suggest widespread reorganization of the 
oceanic and atmospheric circulation during the MWP- 
LIA transition,466-68 leaving the question open of 
the climate state, forcing, and teleconnections dur- 
ing Medieval times. The origin of this reorganization, 
possibly resulting from random internal variations 
or externally driven forcing? or a combination 
thereof, requires further exploration, but will ulti- 
mately contribute to a dynamical understanding of 
paleoclimatic variability and its forcing. 

With an ultimate goal of understanding the 
climate system well enough to make predictions, com- 
parisons of reconstructed and modeled climate varia- 
tions are necessary. To provide a true test of skill, these 
two sources of information should ideally be kept 
independent until they are both reasonably accom- 
plished. Intermingling at an earlier stage may lead to 
overfitting and thus overly confident assessments of 
predictive capability. It is not clear to us if it is ‘too 
early’, but the fact is that such comparisons have been 
made and show that post-1200 a.D. climate models 
and reconstructions provide roughly compatible pic- 
tures of past temperature variations. This agreement 
suggests that either (1) to the first order the climate 
of the past ~800 years is understood at the hemi- 
spheric scale (2) climate models (or reconstructions) 
have been tweaked until they fit, or, but not necessar- 
ily unrelated to the above, (3) a good fit is obtained 
for the wrong reasons (e.g., an underestimate of the 
natural forcing and an unrealistically large climate 
sensitivity of the models, or vice versa, may compen- 
sate for each other), or (4) some combination these 
factors. Resolution of this predicament may require a 
bottom-up and community wide re-evaluation of the 
basic data, assumptions, and algorithms. 


WHERE NEXT? 


With respect to the reconstruction of large-scale 
temperature variation, increases in the amount and 
quality of proxy data and improved understanding of 
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how to best combine these data have led to significant 
progress being made in the past 30 years; however, 
by far, not all is solved. Here, we outline a couple of 
additional research strategies that may contribute to 
reducing uncertainties in future reconstructions and 
contribute to improved predictions. 


Expert Assessment 

The diverse nature of proxy climate archives, rang- 
ing from tree-ring isotopes to documentary data to 
marine sediment particle sizes, means that a wide vari- 
ety of physical, biological, and chemical processes and 
their assessment are responsible for the encoding and 
decoding of the climate signal. A detailed familiarity, 
and ideally a mechanistic understanding, of these pro- 
cesses and the characteristics (e.g., seasonal response, 
nonlinearities, anthropogenic disturbances, frequency 
biases) of different archives in general and individ- 
ual records, in particular, is likely advantageous if 
not invaluable for developing climate reconstructions. 
Every reconstruction—ranging from the very simple 
normalization and averaging of all available records 
to those so complicated that GCM models need to be 
invoked to even get an idea if they may work with 
real world data’'—depend fundamentally upon the 
individual data archives. While climate reconstruction 
may be more complicated than problems that could 
be posed to the first mechanical computer by its inven- 
tor, a quote from the 1864 autobiography of Charles 
Babbage,” may still be apt: 


On two occasions I have been asked,—‘Pray, 
Mr. Babbage, if you put into the machine wrong 
figures, will the right answers come out?” In one case 
a member of the Upper, and in the other a member 
of the Lower, House put this question. I am not able 
rightly to apprehend the kind of confusion of ideas 
that could provoke such a question. 


Expert assessment to evaluate the signal of a particular 
record from a particular proxy archive (e.g., the low- 
frequency skill of a new speleothem record) will 
be invaluable in trying to minimize ‘wrong figures’ 
being put into a large-scale reconstruction. It seems 
advisable at this point to use fewer, but expert- 
assessed proxy records, rather than hundreds of proxy 
series, and hope that reconstruction algorithms will 
overcome the often huge noise components typical 
for many of the available time series. This seems 
particularly to be the case for the lowest frequency 
variations and the lower resolution proxies that 
cannot even be calibrated (and verified) against the 
short instrumental data. Specific considerations from 
experts acquainted with the proxy archives will likely 
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vary considerably from field to field and also as 
advances are made in quantifying and understanding 
proxy signals, and as new proxy types, parameters, 
and processing methods are developed. 


Ensemble Reconstructions 

Although high-quality data are keys, new methods 
that explicitly recognize and incorporate uncertain- 
ties inherent to the reconstruction process could 
lead to the development and implementation of 
ensemble/probabilistic assessments of past climatic 
variation? (But see also Refs 20,73,74). It is unlikely 
that structural uncertainties, such as the use of a 
single reconstruction methodology, a particular cal- 
ibration model, a preferred approach to preserve 
low-frequency variations in tree-ring data, and selec- 
tion of some truncation parameter in multivariate 
regularization approaches, will be eliminated. Intu- 
ition for the characteristics of the underlying problem 
may limit this range of choices to a reasonable, but not 
unique set of possibilities. Systematic exploration and 
acceptance of these ensemble members as viable alter- 
natives may minimize any subconscious selection of 
outcomes that fit with preconceived notions”? for how 
the climate has varied. A simple example is for users 
to simultaneously consider the suite of reconstruc- 
tions shown in Figure 1, rather than selecting and only 
reporting results based upon agreement with a partic- 
ular reconstruction. In practice, more comprehensive 
schemes can be developed that consider uncertainty 
and decisions at numerous levels (e.g., data uncer- 
tainties, pre-processing steps, data selection, method 
definition). However, the ensemble range should avoid 
useless outcomes and thus be small enough to convey 
discriminant power. This is one reason why a mech- 
anistic understanding and expert assessment of the 
available proxies is invaluable. 


Forcing Time-series 

When linking reconstructed and modeled past temper- 
atures to estimate a value for the climate sensitivity or 
the efficacy of the more poorly quantified past forcings 
(namely, volcanic,”° solar,” and land-use/land-cover 
changes’*) more than accurate temperature recon- 
structions are necessary.” This is readily apparent 
when it is considered that the multidecadal surface 
temperature change averaged over the NH in even the 
most complex GCMs can largely be boiled down to 
a linear transformation of the input forcing. Climatic 
effects of individual volcanic eruptions are usually 
short-lived, but clustering of eruptions, as (inconve- 
niently for forcing discrimination) occurred during 
the Late Maunder and Dalton solar minima, can also 
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induce apparent low-frequency climatic variations.®° 


Land-use shifts may elicit strong climate changes, 
but this response is itself unclear, and climate mod- 
els simulations driven by the same land-use history 
display a large range of responses.*! In a similar 
vein, state of the art simulations®* driven by cur- 
rent best estimates (reduced by a factor of 2 from 
the 2001-2007 IPCC reports) of solar forcing fail to 
reproduce the MWP-LIA sequence prominent in most 
proxy records,” and therefore may require explanation 
from internal random unforced variations, attribution 
to volcanic activity, or other poorly parameterized 
factors. Control simulations with constant external 
forcing produce multidecadal variations of the NH 
temperature of just 0.2°C,*° seemingly ruling out sole 
effects of internal random climate variations. Given 
the defined relationships between forcing series and 
modeled output, refined forcing series®?:’®784 may 
be crucial to resolve such questions. Constraining the 
large-scale amplitude and efficacy of forcing agents 
appears to be the highest priority, although as models 
are increasingly being used for regional studies, par- 
allel increases in spatial and temporal resolution of 
forcing time-series will also be necessary. 


CONCLUSION 


The geological principle of uniformitarianism, often 
summarized as ‘the present is the key to the past’, 
was invoked by early natural scientists to comprehend 
long time-scale processes that occurred so slowly or 
in the far distant past, that only traces remaining 
could be used to infer what occurred. Significant 
natural variability and the long response time of the 
oceanic and terrestrial systems to climate change, 
make it exceedingly difficult to understand and 
therefore predict climate based upon short present 
observational data. Current anthropogenic activities 
have led to unprecedented trajectories and states in the 
earth’s coupled climate system, but characterization 
of the natural climate variability will at least allow a 
better understanding of the basic operating rules and 
patterns of climate change. In the context of trying to 
understand the consequences of new anthropogenic 
regimes, this uniformitarianism paradigm must be 
flipped and the past used to grasp hints for what 
the future may hold. 

We have herein reviewed the changing icons 
in high-resolution paleoclimateology and discussed 
the current (post-2007 IPCC) understanding of large- 
scale temperature variation over the past 1000 or 
so years. We suggest that expert knowledge of the 
proxy archives will be critical to reduce uncertainties 
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and large-scale reconstructions. At the same time, 
recognition and acceptance of structural uncertainties 
will allow a more accurate assessment of reconstruc- 
tion errors. Improved forcing time series will allow 
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